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cn ■ ABSTRACT 
\D '. 

O ' Main sequence turnoff magnitudes from the recent set of Yale isochrones (Chaboyer 

g : et al. 1995) have been combined with a variety of relations for the absolute magnitude 

\ of RR Lyr stars (Mv(RR)) to calibrate age as a function of the difference in magnitude 

^ ■ between the main sequence turn-off and the horizontal branch (AV^b)- A best esti- 

^ . mate for the calibration of Mv(RR) is derived from a survey of the current literature: 

Mv(RR) = 0.20 [Fe/H] + 0.98. This estimate, together with other calibrations (with 
P \ slopes ranging from 0.15 to 0.30) has been used to derive AVj^b ages for 43 Galactic 

globular clusters. Independent of the choice of Mv(RR), there is no strong evidence for 
C3 ■ an age-Galactocentric distance relationship among the 43 globular clusters. However, an 

age-metallicity relation exists, with the metal-poor clusters being the oldest. A study 
^ of the age distribution reveals that an age range of 5 Gyr exists among the bulk of 

■ the globular clusters. In addition, about 10% of the sample are substantially younger, 

^ and including them in the analysis increases the age range to 9 Gyr. Once again, these 

statements are independent of the Mv(RR) relation. Evidence for age being the second 
parameter governing horizontal branch morphology is found by comparing the average 
AV^B age of the second parameter clusters to the normal clusters. The second parame- 
ter clusters are found to be on average 2-3 Gyr younger than the other clusters, which 
is consistent with age being the second parameter. These results suggest that globular 
clusters were formed over an extended period of time, with progressively more metal-rich 
globular clusters ([Fe/H] ;^ —1-7) being formed at later times. 

Subject headings: globular clusters: general - Galaxy: formation - Galaxy: Halo 



^NOAO is operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science 
Foundation. 
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1. Introduction 

Understanding tiie process of galaxy formation 
continues to be one of the key quests in astrophysics. 
In this regard, the Milky Way plays a unique role be- 
cause it is the only galaxy for which we can obtain 
detailed chemical, kinematic and chronological infor- 
mation. Observational and theoretical studies over 
the last 60 years have lead to a basic understand- 
ing of how the Galaxy formed (see Larson 1991). It is 
clear that the spherical, metal-poor halo of our galaxy 
formed early during the collapse of the proto-galactic 
cloud. The collision at the mid-plane halted the gas 
collapse and lead to the formation of the rotating, 
thin disk. 

However, there are many unanswered questions re- 
garding the formation of the Galaxy. When did the 
bulge form? How and when did the thick disk form? 
How important is later infall and accretion? Did the 
halo form over an extended period of time? Was halo 
formation a chaotic or smooth process? An important 
step towards answering these questions is to deter- 
mine accurate ages for the various stellar populations. 
Globular clusters (GCs) play a key role in this regard, 
for their derived ages are the most accurate of any ob- 
ject in the halo and thick disk/bulge. In this paper, 
ages for 43 Galactic GCs which have well observed 
colour magnitude diagrams (CMDs) are derived and 
analyzed to probe the formation of the Galactic halo. 

Information regarding the formation and evolution 
of the Galaxy has traditionally been obtained by sur- 
veys of stars or star clusters which have one or more of 
the following properties measured: locations, metal- 
licities, velocities, and ages. The classic paper by 
Eggen, Lynden-Bell & Sandage (1962) analyzed ul- 
traviolet excesses, radial velocities and proper mo- 
tions of nearby stars to conclude that the Galactic 
halo formed during the rapid, monolithic collapse of 
the proto- Galactic gas cloud. Evidence for a quite dif- 
ferent halo formation theory was presented by Searle 
& Zinn (1978). On the basis of their studies of GC 
metallicities and horizontal branch morphology Searle 
& Zinn (1978) proposed that the halo formed via ac- 
cretion over several gigayears (Gyr) in a rather chaotic 
manner. These contrasting theories continue to be a 
critical part of discussions of Galactic formation (see 
Majewski 1993 for a recent review). In this regard, 
the determination of accurate absolute and relative 
ages for GCs plays an important role in discovering 
the time scale for the formation of the Galactic halo. 

There are a variety of different methods which can 
be used to derive the ages of GCs. All of these tech- 
niques rely on comparing some aspect of an observed 



CMD to theoretical stellar models or isochrones. The 
most accurate relative ages can be derived using the 
difference in colour between the main sequence turn- 
off and the base of the red giant branch (A(B — V), 
Sarajedini & Demarque 1990; VandcnBerg. Bolte & 
Stetson 1990). However, the colours of theoretical 
isochrones are very uncertain, as they depend on stel- 
lar atmospheres and the mixing length treatment of 
convection. As such, transforming an observed dif- 
ference in A(B — V) into an age difference is siibject 
to large theoretical uncertainties. These uncertainties 
can be minimized by comparing clusters with similar 
metallicities. 

If one wishes to inter-compare ages of clusters with 

different metallicities, then the difference in magni- 
tude between the main sequence turn-off and the hor- 
izontal branch (AV^b) yields ages which have the 
smallest theoretical errors. Unfortunately, it is of- 
ten difficult to determine AV^g observationally, so 
that the error in the derived age can be rather large 
(^ 10%). The absolute magnitude of the main se- 
quence turn-off (Mv(TO)) is a well determined the- 
oretical quantity. The new set of Yale isochrones 
(Chaboyer et al. 1995) provide an up-to-date calibra- 
tion of My (TO) for a wide range of ages and chemical 
compositions. The absolute magnitude of the hori- 
zontal branch (HB) is independent of age (over the 
range ~ 8 ~ 22 Gyr), however its absolute level is 
not well determined in theoretical models due to the 
importance of convection and semi-convection in the 
nuclear burning regions of these stars. Fortunately, 
there are a variety of independent, observationally 
based methods which can be used to determine the 
absolute magnitude of RR Lyr stars (Mv(RR)) which 
lie on the HB. Hence, the AV^g ages derived in this 
paper are based on the calibration of Mv(TO) 
function of age and metallicity from the new set of 
Yale isochrones, coupled with a variety of determina- 
tions of Mv(RR). 

There have been a number of studies of AVyg ages 
of GCs in recent years (Sarajedini & King 1989; 
Sandage & Cacciari 1990; Carney, Storm & Jones 
1992; Chaboyer, Sarajedini & Demarque 1992; Walker 
1992a; Caputo et al. 1993; Sandage 1993). These 
studies have shown that the choice of a Mv(RR) re- 
lation is crucial to the conclusions which are drawn 
based on AV^^ ages. This work differs from previ- 
ous studies in three main ways: (i) use of the new 
Yale isochrones to determine Mv(TO); (ii) an ex- 
panded observational database of observed AV^g val- 
ues (30% more than in our 1992 compilation); and 
(iii) the use of a large number of Mv(RR) relations 
to explore in detail how the choice of Mv(RR) affects 
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our conclusions. 

A brief description of the new Yale isochrones along 
with a discussion of Mv (RR) and the theoretical cal- 
ibration of AV^g is presented in §2. Section 3 re- 
views the basic observational data and tabulates the 
AVyg ages. A discussion of the correlations between 
age, metallicity and galactocentric distance is con- 
tained in §4. Evidence for an age range within the 
Galactic globular cluster system is presented in §5. 
Section 6 examines the second parameter problem in 
the context of the AV^g ages. Finally, §7 discusses 
the major results of this paper, and their implications 
for the formation of the Galactic halo. 

2. Theoretical Calibration of AVj^g 
2.1. Mv(TO) 

The recent set of Yale isochrones (Chaboyer et al. 
1995) are used to provide a calibration of Mv(TO) as 
a function of age and metal abundance. These iso- 
chrones are based on new stellar evolution models 
which incorporate the latest available input physics: 
opacities from Iglesias & Rogers (1991, high tem- 
perature) and Kurucz (1991, low temperature) and 
nuclear reaction rates from Bahcall & Pinsonneault 
(1992) and Bahcall (1989). The colour transforma- 
tion of Green, Demarque & King (1987) was used to 
construct isochrones in the observational plane. The 
new set of Yale isochrones are based on standard stel- 
lar models, which do not include the effects of diffu- 
sion, or the Debye-Hiickel correction to the equation 
of state. Including these two effects would system- 
atically reduce the GC ages presented in Table 3 by 
~ 13% (Chaboyer 1995). 

In order to span the range of metallicities of ob- 
served globular clusters, Mv(TO) values were deter- 
mined from isochrones with [Fe/H] = —2.8,-2.3, 
-1.8, -1.3, -1.0 and -0.44. The isochrones with 
[Fe/H] < —1.0 have a helium abundance oiY = 0.23. 
This is in good agreement with recent determinations 
of the primordial helium abundance (Pagel & Ka- 
zlauskas 1992; Balbes, Boyd & Mathews 1993; Izo- 
tov, Thuan & Lipovetsky 1994). The most metal 
rich isochronc ([Fc/H] — —0.44) has a helium abun- 
dance oi Y = 0.25 and assumes that the a-capture 
elements (O, Mg, Si, S, and Ca) are enhanced by 0.2 
dex over their solar values (ie. [a/Fe] — -1-0.2). The 
more metal-poor isochrones assume [a/Fe] = -1-0.4. 
These [a/Fe] values were chosen to be in agreement 
with observations of halo stars (Lambert 1989; Dick- 
ens et al. 1991; King 1994; Nissen et al. 1994). The 
new Yale isochrones are tabulated every 2 Gyr for the 
older ages (10 - 22 Gyr). In order to provide a finer 



grid for interpolation purposes in this study, we have 
recomputed these isochrones using a 1 Gyr spacing, 
between 8-22 Gyr. 

2.2. Mv(RR) 

There are numerous observational and theoretical 
techniques which may be used to derive Mv(RR). 
There is general agreement that the absolute mag- 
nitude of RR Lyr stars is given by an equation of the 
form 

Mv(RR) = Ai[Fe/H] -f 7, (1) 

where /i is the slope with metallicity and 7 is the 
zero-point. The zero-point is important for setting 
the overall absolute ages, while the slope is impor- 
tant in determining the relative ages for GCs with 
different metallicities. Some techniques for deter- 
mining Mv(RR) are best for determining the zero- 
point, while other techniques are best at deriving the 
slope. Due to possible systematic efi'ects, the Baade- 
Wesselink and infrared flux analysis are best used to 
determine the slope with metallicity. Such analyses 
of field RR Lyr stars have been published recently by 
two groups. Jones et al. (1992) found = 0.16±0.03, 
while Skillen et al. (1993) determined /.t = 0.21±0.05. 
The theoretical HB models of Lee (1990) should also 
give a reliable determination of the slope, and yield 
/I — 0.18 ± 0.01. From an analysis of the Ooster- 
hoff period shift effect in GCs, Sandage (1993)[| found 
^jL = 0.30 ± 0.12. Thus, it appears that Mv(RR) has 
a rather shallow slope with metallicity of ^ ~ 0.20, 
though it may be somewhat premature to totally ex- 
clude slopes as high as ^ = 0.30. 

A reliable determination of the zero-point in equa- 
tion (^ can be made by measuring the apparent mag- 
nitude of a number of RR Lyr stars in the LMC and 
then using the distance to the LMC to obtain 7. This 
is the approach used by Walker (1992a) who found 
Mv(RR) = 0.44 ±0.10 at [Fe/H] = -1.9 (implying 
7 = 0.82±0.10), assuming (m - M)lmc = 18.5±0.10. 
This distance modulus to the LMC was based on main 
sequence fitting, and analysis of the Cepheid variables 
and the rings associated with SN1987A. However, 
Gould (1995) has recently re-analyzed the SN1987A 
distance estimate, and determined an upper limit of 
(m — M)lmc = 18.37. Using this distance estimate 
to the LMC and Walker's (1992a) RR Lyr photom- 
etry, one finds Mv(RR) = 0.57 at [Fe/H] = -1.9 

^Sandage (1993) gives an error of 0.12 for the slope of Mboi (RR) 
with metallicity, but does not quote an error in Mv (RR) . How- 
ever, Sandage obtains Mv(RR) from a simple transformation 
of his Mboi (RR) equation (Mboi = Mv -I- 0.06 [Fe/H] -|- 0.06), 
so we have taken the error in his Mv(RR) relation from his 
quoted error in his M^oi (RR) relation. 
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Table 1 
Fitting Coefficients 



Mv(RR) 








Standard 










Red HB 






/3o 


Pi 








ft 


/3o 


/3i 




ft 


ft 




0.17 [Fe/H] + 


0.79 


74.995 


-47.943 


8.435 


6.219 


0.551 


-2.025 


100.913 


-62.728 


10.550 


21.501 


3.562 


-5.436 


0.20 [Fe/H] + 


0.98 


77.703 


-51.506 


9.419 


4.669 


0.457 


-1.611 


88.669 


-58.294 


10.485 


18.692 


3.385 


-4.810 


0.15 [Fe/H] + 


0.98 


83.026 


-54.407 


9.813 


7.659 


0.573 


-2.664 


81.215 


-53.656 


9.768 


22.377 


3.664 


-6.078 


0.15 [Fe/H] + 


0.72 


76.664 


-48.252 


8.335 


7.177 


0.614 


-2.297 


105.035 


-64.100 


10.550 


23.136 


3.675 


-5.858 


0.20 [Fe/H] + 


1.06 


78.270 


-52.404 


9.725 


5.012 


0.456 


-1.752 


69.791 


-47.837 


9.144 


19.122 


3.474 


-5.006 


0.20 [Fe/H] + 


0.82 


72.434 


-46.790 


8.343 


4.684 


0.503 


-1.499 


99.041 


-62.095 


10.550 


19.475 


3.408 


-4.803 


0.25 [Fe/H] + 


1.14 


74.978 


-51.177 


9.753 


2.728 


0.406 


-0.880 


85.981 


-57.814 


10.783 


19.077 


3.341 


-4.855 


0.25 [Fe/H] + 


0.91 


67.189 


-44.791 


8.296 


2.212 


0.434 


-0.656 


93.237 


-60.091 


10.550 


16.014 


3.195 


-3.748 


0.30 [Fe/H] + 


1.22 


72.167 


-50.378 


9.842 


0.000 


0.240 


0.000 


81.425 


-56.089 


10.783 


15.884 


3.125 


-3.776 


0.30 [Fe/H] + 


1.01 


67.561 


-45.831 


8.638 


0.000 


0.249 


0.000 


86.174 


-57.216 


10.420 


12.826 


3.024 


-2.718 



(7 = 0.95). Using a statistical parallax analysis 
of field RR Lyr stars, Layden, Hanson & Hawley 
(1994) found a zero-point which is 0.24 mag fainter 
(7 = 1.06) than that quoted by Walker (1992a) (and 
0.11 mag fainter than the revised Walker value above). 
This suggests that 7 ~ 0.98 is a reasonable choice, 
and so our best estimate for the absolute magnitude 
of the RR Lyr stars is Mv(RR) = 0.20[Fe/H] + 0.98. 
This choice of Mv(RR) predicts that the RR Lyrs 
in M92 (with [Fe/H] = -2.24 ± 0.08) should have 
Mv(RR)m92 = 0.53 ± 0.02. In their Baade-Wesselink 
analysis of M92 RR Lyrs, Storm, Carney & Latham 
(1994) determined ]V[v(RR)m92 = 0-43 ± 0.22, which 
is in reasonable agreement. However, the possible 
0.1 mag discrepancy would reduce the ages for the 
metal-poor clusters by 10% and so a more precise de- 
termination of Mv (RR) in metal-poor GCs is clearly 
desirable. Although we have derived our best esti- 
mate for the absolute magnitude of the RR Lyr stars 
(Mv(RR) = 0.20[Fe/H] + 0.98), in order to explore 
in a systematic way the effect that uncertainties in 
Mv(RR) have on the GC ages, ages will be derived 
using both the Walker (1992a) and Layden et al. 
(1994) zero-points with slopes which vary from 0.15 
to 0.30. 

2.3. Derivation of AV^g Ages 

The Mv(TO) values from the new Yale isochrones 
are combined with a given Mv(RR) relation to form 
a grid, which specifies age given AV^g and [Fe/H]. 
The grid is then fit to an equation of the form 

t9 = (3o + 0iAY + (32AY^ + (33[Fe/U] 

+ /34[Fe/H]'+/35AV[Fc/H], (2) 



where tg is the age in Gyr. The rms residuals of the 
points from the fit were about 0.15 Gyr. The above 
formula is used to determine ages for GC which have 
RR Lyrs, or a blue HB. In the case clusters with 
purely blue HBs (i.e. few or no RR Lyrae variables), 
observers quote the V mag of the blue edge of the in- 
stability strip, which is usually a reasonably accurate 
measurement of Mv(RR). In some cases, they com- 
pare to the blue HB of a cluster with RR Lyraes to 
infer Mv(RR). 

In the case of clusters with red HBs (HB typeQ 
< — O.8.), the situation is slightly more complicated. 
In these clusters, observers usually quote the mean 
or median mag of the red HB stars. This quantity 
can be anywhere from 0.05 to 0.2 mag brighter or 
fainter than the RR Lyr level depending on the clus- 
ter metallicity and age. In order to correct for this 
effect, a semi-empirical approach is taken. The off- 
set between the red HB level and RR Lyr level may 
be determined from theoretical HB models, and this 
correction can then be applied to the red HB clusters, 
as discussed by FuUton et al. (1995). As this offset 
depends on relative quantities in the theoretical mod- 
els, it should be reasonably reliable. HB models by 
Lee, Demarque & Zinn (1987), Dorman (1992) and 
Castellani, Chieffi & Pulone (1991) find offsets which 
agree to within 0.05 mag. The offsets used in this 
study are derived from HB models kindly provided 
to us by Lee (1995). Lee constructed synthetic HB 
models of red HB clusters with a range of ages and 
abundance from which he has calculated Mv(HB). 

^The HB type has the following definition: HBtype = (B — 
R)/(B -I- V + R), where B is the number of HB stars blueward 
of the instability strip, R is the number of HB stars redward of 
the instability strip, and V is the number of RR Lyr stars 
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Given the relation Mv(RR) 0.17 [Fe/H] + 0.79, 
which comes from the Lee HB models, the function 
5 = Mv(redHB) — Mv(RR) can be calculated as a 
function of [Fe/H] and age. For other RR Lyrae lu- 
minosity relations, the function S is used to correct 
Mv(RR) to Mv(redHB). 

The coefficients of the fit (equation ^ for the 
standard and red HB cases are given in Table 1. 
These coefficients have been tabulated for 10 dif- 
ferent Mv(RR) relations: the Lee (1995) relation 
(Mv(RR) = 0.17 [Fe/H] +0.79); our best estimate 
for the true relation (Mv(RR) = 0.20 [Fe/H] + 0.98); 
and relations with slopes of 0.15, 0.20, 0.25 and 0.30 
using the zero-points of Walker (1992a) and Layden 
et al. (1994). These Mv(RR) relations span the 
range reported by various groups using a variety of 
observational and theoretical techniques (see §2.2). 
The coefficients presented in Table 1 are valid for 
—2.8 < [Fe/H] < —0.44 and for ages in the range 
8-22 Gyr. 

3. The Ages 

Estimating GC ages using AV^g requires an accu- 
rate measurement of V(TO) and V(HB), along with 
an estimate of [Fe/H]. Continued advances in CCD 
technology and image reduction, along with the ad- 
vent of the HST have lead to a wealth of high qual- 
ity CMDs of GCs in recent years. Table 2 lists var- 
ious observational quantities for 43 GCs for which 
reliable age determinations may be made based on 
published observations of V(TO) and V(HB). To be 
conservative, M22 and w Cen have not been included 
in this group, as there is evidence for a range in 
metallicity in these clusters (Noble et al. 1991; Smith 
1987) which complicates the age determination pro- 
cess. References for V(TO) and V(HB) are provided 
in the table. In some cases, the observers do not 
quote V(HB); rather, they provide the apparent mag- 
nitude of the zero-age horizontal branch (V(ZAHB)). 
These have been converted to mean HB magnitudes 
using equation 4 from Carney et al. (1992: V(HB) = 
V(ZAHB) - 0.05 [Fe/H] - 0.20). The V(HB) (and cor- 
responding AVyg values) which have been corrected 
for this effect are indicated by an asterisk next to 
the V(HB) value in Table 2. Our 1992 compilation 
(Chaboyer et al. 1992) does NOT include a correc- 
tion for this effect, as we were not aware that some 
observers quoted V(ZAHB). 

The cluster [Fe/H] values and their errors are taken 
from Zinn & West (1984), except for NGC 4147, 5053 
(Armandroff, Da Costa & Zinn 1992), NGC 6218, 
Ter 7, Ter 8, and Arp 2 (Da Costa & Armandroff 



1995), and Rup 106 (Da Costa, Armandroff, & Norris 
1992). The reddennings are from Zinn (1985), except 
for NGC 4590 (Walker 1994), NGC 5053 (Sarajedini 
& Milone 1995), NGC 6352 (FuUton et al. 1995), 
NGC 6535 (Sarajedini 1994), NGC 6584 (Sarajedini 
& Forrester 1995), NGC 6652 (Ortolani, Bica & Bar- 
buy 1994), Ter 7 (Webbink 1985), Ter 8 (Ortolani 
& Gratton 1990), Arp 2 (Buonanno et al. (1995a) 
and Rup 106 (Da Costa, Armandroff & Norris 1992). 
The Galactic coordinates for the clusters are taken 
from Shawl & White (1986), except for Rup 106 and 
Pal 12, which are from Webbink (1985). The HB 
types are from Lee, Demarque & Zinn (1994), with 
the exceptions of NGC 4590 (Walker 1994), NGC 
6584 (Sarajedini & Forrester 1995), NGC 6535 (Sara- 
jedini 1994), NGC 6652 (Ortolani et al. 1994), Ter 
7, Arp 2 (Buonanno et al. 1995a,b), Ter 8 (Or- 
tolani & Gratton 1990), and IC4499 (Ferraro et al. 
1995). The groupings into disk, old halo and younger 
halo clusters are from Zinn & Lee (1995; see also 
Zinn 1993). In calculating the Galactocentric distance 
(Rgc) of each cluster, we have adopted Rq = 8.0 kpc, 
Ay — 3.2E(B — V), and a distance modulus de- 
rived from V(HB) and our preferred Mv(RR)relation 
(Mv(RR) = 0.20 [Fe/H] -I- 0.98). Proper motion stud- 
ies exist for 16 of the GCs in our sample, and these 
have been used by Majewski (1994) to determine the 
apogalactica distances (Rapo) listed in Table 2. As ex- 
pected, most of the apogalactica distances are quite 
similar to the Galacto-centric distances. However, 
there are a few notable exceptions. The present po- 
sitions of NGC 5466, 5904 and 6205 are considerably 
smaller than their apogalactica distances. 

Using the AV^g and [Fe/H] values listed in Ta- 
ble 2, GC ages are derived using equation (^), with 
the (3 coefficients listed in Table 1. The error in the 
derived age is calculated by propagating the errors 
in AVj![B ^"^^ [Fe/H] through equation (||). The er- 
ror in the derived age is dominated by the error in 
AVyg. For the statistical analysis which comprises 
the bulk of this paper, it is important that the error in 
AVyg represent a Gaussian 1-sigma error bar. How- 
ever, it is doubtful that the observers quote such an 
error bar. In order to get an estimate for the Gaussian 
1-sigma error bar, the literature has been searched 
for independent measurements of AVj^g. Appendix 
A presents an analysis of these independent observa- 
tions, and concludes that a reliable estimate for the 
Gaussian 1-sigma error in AV^g may be obtained by 
multiplying the quoted error by 0.61. This correction 
factor has been applied to all of the quoted AV^g er- 
rors when determining the error in the derived age. 
In some cases, errors in AV^g were not given by the 
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Table 2 

Globular Cluster Parameters 



Cluster 


e b 


Rgc 


R-apo 


References'^ 


NGC Name 


(degrees) (degrees) E(B - V) [Fc/H] 


V(HB)=» AVTOb (kpc) 


(kpc) 


HB type Group<= HB TO 



104 


47 Tug 


305.896 


-44.900 


0.04 


-0.71 ± 0.08 


14, 


288 




151.328 


-89.383 


0.04 


-1.40 ± 0.12 


15, 


362 




301.533 


-46.248 


0.06 


-1.27 ± 0.07 


15, 


1261 




270.541 


-52.126 


0.00 


-1.31 ± 0.09 


16, 


1851 




244.514 


-35.037 


0.02 


-1.36 ± 0.09 


16, 


1904 


M79 


227.229 


-29.351 


0.01 


-1.69 ± 0.09 


16, 


2298 




245.628 


-16.007 


0.08 


-1.85 ±0.11 


16, 


2808 




282.192 


—11.253 


0.22 


-1.37 ±0.09 


15. 


3201 




277.228 


8.641 


0.21 


—1.61 ± 0.12 


14. 


4147 




252.850 


77.189 


0.02 


-1.80 ±0.14 


17, 


4590 


M68 


299.626 


36.052 


0.07 


—2.09 ± 0.11 


15, 


5024 


M53 


332.965 


79.765 


0.00 


-2.04 ±0.08 


16, 


5053 




335.695 


78.944 


0.06 


-2.41 ±0.06 


16, 


5272 


M3 


42.208 


78.708 


0.01 


-1.66 ±0.06 


15, 


5466 




42.150 


73.592 


0.00 


-2.22 ± 0.36 


16, 


5897 




342.946 


30.294 


0.06 


-1.68 ± 0.11 


16, 


5904 


M5 


3.863 


46.796 


0.03 


-1.40 ± 0.06 


14, 


6101 




317.747 


-15.825 


0.04 


-1.81 ± 0.15 


16, 


6121 


M4 


350.974 


15.972 


0.40 


-1.33 ± 0.10 


13, 


6171 


M107 


3.374 


23.011 


0.31 


-0.99 ± 0.06 


15, 


6205 


M13 


59.008 


40.912 


0.02 


-1.65 ± 0.06 


14, 


6218 


M12 


15.717 


26.313 


0.17 


-1.34 ± 0.09 


14, 


6254 


MIO 


15.138 


23.077 


0.32 


-1.60 ± 0.08 


14, 


6341 


M92 


68.339 


34.860 


0.02 


-2.24 ± 0.08 


14, 


6352 




341.421 


-7.167 


0.21 


-0.51 ± 0.08 


15, 


6397 




338.165 


-11.958 


0.18 


-1.91 ± 0.14 


12, 


6535 




27.177 


10.436 


0.44 


-1.75 ±0.15 


15, 


6584 




342.144 


-16.414 


0.07 


-1.54 ±0.15 


16, 


6652 




1.534 


—11.377 


0.10 


—0.89 ± 0.15 


15, 


6752 




336.496 


-25.627 


0.04 


-1.54 ±0.09 


13, 


6809 


M55 


8.795 


-23.270 


0.06 


-1.82 ±0.15 


14, 


6838 


M71 


56.744 


-4.564 


0.27 


-0.58 ± 0.08 


14, 


7006 




63.770 


-19.407 


0.05 


-1.59 ±0.07 


18, 


7078 


M15 


65.013 


-27.313 


0.10 


-2.15 ±0.08 


15, 


7099 


M30 


27.180 


-46.835 


0.04 


-2.13 ±0.13 


15, 


7492 




53.386 


-63.478 


0.00 


-1.82 ±0.30 


17, 




Ter 7 


3.387 


-20.063 


0.06 


-0.36 ±0.09 


17, 




Ter 8 


5.758 


-24.558 


0.20 


-1.99 ±0.08 


17, 




Rup 106 


300.888 


11.670 


0.20 


-1.69 ±0.05 


17, 




Pal 5 


0.852 


45.860 


0.03 


-1.47 ±0.29 


17, 




Pal 12 


30.510 


-47.680 


0.02 


-1.14 ±0.20 


17, 




IC4499 


307.354 


-20.473 


0.25 


-1.75 ±0.20 


17, 




Arp 2 


8.543 


-20.787 


0.08 


-1.70 ±0.11 


18, 



,09 


3.61 ± 0.10 


7.3 


7.3 


-1.00 ± 0.03 


D 


1 


1 


,27* 


3.73 ± 0.12 


11.2 


11.5 


0.95 ± 0.08 


OH 


2 


2 


,29' 


3.56 ± 0.14 


8.8 


9.2 


-0.87 ±0.08 


YH 


2 


2 


,57* 


3.57 ± 0.12 


16.8 




-0.70 ± 0.10 


YH 


3 


3 


,15 


3.45 ± 0.10 


16.5 


— 


-0.33 ± 0.08 


YH 


4 


4 


,03* 


3.57± — 


17.8 




0.89 ± 0.16 


OH 


5 


5 


,11 


3.49 ±0.21 


16.1 


— 


0.93 ± 0.24 


OH 


6 


6 


,97* 


3.63 ±0.14 


10.2 




-0.54 ± 0.06 


YH 


2 


2 


,76 


3.39 ±0.17 


8.8 


— 


0.08 ± 0.06 


YH 


7 


8 


,00 


3.60 ± — 


20.5 


25.7 


0.55 ±0.14 


YH 


9 


9 


,64 


3.42 ±0.10 


9.6 




0.44 ± 0.05 


YH 


10 


10 


,90 


3.55 ±0.14 


18.9 


— 


0.76 ±0.10 


OH 


11 


12 


,65 


3.48 ± — 


16.2 


— 


0.61 ±0.18 


OH 


13 


14 


,63 


3.54 ± 0.09 


11.7 


18.3 


0.08 ±0.04 


YH 


15 


15 


,62 


3.58 ± — 


16.7 


29.2 


0.68 ± 0.14 


OH 


16 


17 


,35 


3.60 ±0.18 


7.6 


— 


0.91 ± 0.10 


OH 


18 


18 


,98* 


3.62 ± 0.11 


6.0 


21.1 


0.37 ± 0.06 


OH 


2 


2 


,60 


3.40 ± — 


10.7 


— 


0.84 ± 0.16 


OH 


19 


19 


,22* 


3.68 ± 0.16 


6.4 


6.4 


-0.07 ±0.10 


OH 


2 


2 


,55* 


3.75 ± 0.18 


3.6 


3.9 


-0.76 ± 0.08 


OH 


2 


2 


,83* 


3.67 ±0.21 


8.2 


18.9 


0.97 ± 0.08 


OH 


2 


2 


,90 


3.45 ± — 


4.3 


4.5 


0.92 ± 0.10 


OH 


20 


21 


,65 


3.75 ± 0.15 


4.9 


— 


0.94 ± 0.10 


OH 


22 


23 


,96* 


3.74 ± 0.12 


9.1 


9.3 


0.88 ± 0.08 


OH 


2 


2 


,13 


3.67 ±0.10 


3.6 


— 


-1.00 ± 0.04 


D 


1 


1 


,90* 


3.74 ±0.14 


6.1 


6.1 


0.93 ±0.10 


OH 


2 


2 


,73 


3.66 ±0.19 


4.2 


— 


1.00 ± — 


OH 


24 


24 


,53 


3.47 ± — 


6.9 


— 


-0.09 ± 0.06 


YH 


25 


25 


,85 


3.35 ±0.16 


1.9 




-1.00 ± — 


OH 


26 


26 


,63* 


3.77 ±0.16 


5.4 




1.00 ± 0.04 


OH 


2 


2 


,24* 


3.66 ±0.10 


4.1 




0.91 ±0.10 


OH 


2 


2 


,44 


3.56 ± 0.09 


6.8 


6.7 


-1.00 ±0.04 


D 


27 


27 


,80 


3.55 ±0.12 


36.8 




-0.11 ±0.06 


YH 


28 


28 


,77* 


3.63 ±0.16 


9.9 


10.2 


0.72 ±0.10 


OH 


2 


2 


,11* 


3.62 ±0.14 


6.9 




0.88 ±0.12 


OH 


2 


2 


,52* 


3.72 ±0.14 


23.2 




0.90 ±0.18 


OH 


2 


2 


,76 


3.20 ±0.12 


14.3 




-1.00 ± — 


YH 


29 


29 


,85 


3.65 ± — 


14.3 




1.00 ± — 


OH 


30 


31 


,73* 


3.32 ± 0.07 


17.0 




-0.82 ± 0.15 


YH 


32 


32 


,27* 


3.-53 ±0.14 


15.4 


16.9 


-0.40 ± 0.20 


YH 


2 


2 


,13 


3.30 ± — 


15.1 




-1.00 ±0.12 


YH 


33 


33 


,80 


3.25 ±0.15 


15.7 




0.11 ±0.36 


YH 


34 


34 


,18* 


3.40 ±0.10 


21.5 




0.86 ± — 


OH 


35 


35 



^Entries with an asterick (*) were origininally published as V(ZAHB). The V(HB) and AV^g values have been corrected using V(HB) = 
VzAHB - 0.05[Fe/H] - 0.20 (eqn. 4 from Carney et al. 1992). 

''Errors given in AV^g are those quoted in the original papers, and are not 1 c Gaussian errors. As discussed in Appendix A, a reasonable 
estimate of the 1 a Gaussian error may be obtained by multipling the quoted errors by 0.61 

<=D = Disk; OH = Old Halo; YH = Younger Halo 

''(1) PuUton et al. 1995; (2) Buonanno, Coris & Fusi Pecci 1989; (3) Ferraro et al. 1993; (4) Walker 1992b; (5) Ferraro et al. 1992a; (6) 
Janes & Heasley 1988; (7) Alcaino, Liller & Alvarado 1989; (8) Brewer et al. 1993; (9) Fricl, Heasley & Christian 1987; (10) Walker 1994; 
(11) Cuffcy 1965; (12) Heasley & Christian 1991; (13) Sarajcdini & Milonc 1995; (14) Fahlman, Richer & Ncmcc 1991; (15) Buonanno et al. 
1994; (16) Ncmcc & Harris 1987; (17) Peterson 1986; (18) Sarajcdini 1992; (19) Sarajcdini & Da Costa 1991; (20) Racine 1971; (21) Sato, 
Richer & Fahlman 1989; (22) Hurley, Richer & Fahlman 1989; (23) Harris, Racine & dcRoux 1976; (24) Sarajcdini 1994; (25) Sarajcdini & 
Forrester 1995; (26) Ortolani, Bica & Barbuy 1994; (27) Hodder et al. 1992; (28) Buonanno et al. 1991; (29) Buonanno et al. 1995b; (30) 
Da Costa & Armandroff 1995; (31) Ortolani & Gratton 1990; (32) Buonanno et al. 1993; (33) Stetson et al. 1989; (34) Ferraro et al. (1995); 
(35) Buonanno et al. 1995a 
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observers. For these clusters, an error of 0.083 mag 
was assumed. This is the average of the Gaussian 1- 
sigma errors for those clusters with quoted errors in 

Table 3 presents ages for the 43 GCs, using the 
10 different Mv(RR) relations given in Table 1. The 
heading of each column gives the Mv(RR) relation 
used to derive the ages in that column. These ages 
will be analyzed in detail in the following sections. 
Here, we simply note that ages derived using the Lay- 
den et al. (1994) zero-point for Mv(RR) are approx- 
imately 25% larger than the ages derived using the 
Walker zero-point (e.g. compare columns 5 and 6). 
This illustrates the well known result that a 0.25 mag 
uncertainty in the distance modulus translates into a 
25% uncertainty in ages derived using AV^g. It is 
also interesting to note that there are several young 
clusters in the sample; IC4499, Arp 2, Pal 12, Rup 
106, and Ter 7 have all been shown to be young by 
the A(B — V) technique (see Buonanno et al. 1994) 
and indeed, the AVjSjg ages for these clusters arc all 
young compared to the mean age. NGC 6652, which 
has a small AV^g value (Ortolani, Bica & Barbuy 
1994), also appears to be yoimg. 

Among these young clusters, Ter 7 stands out with 
an age of ~ 9 Gyr, which is at least 2 Gyr younger 
than the others. It appears that Ter 7 is associ- 
ated with the recently discovered Sgr dwarf spheroidal 
galaxy (Ibata, Gilmore & Irwin 1994). Indeed Ter 7, 
Ter 8, Arp 2, and M54 (not on our list) arc located at 
approximately the same distance and the same region 
of the sky as Sgr. In addition. Da Costa & Arman- 
droff (1995) have shown that the above clusters have 
similar radial velocities to that of Sgr. However, an 
inspection of Table 2 reveals that Ter 7 and Ter 8 
have Rgc = 14.3 kpc while Arp 2 has Rgc = 21.5 
kpc. Using the reddening, metallicity and V(HB) 
listed by Da Costa & Armandroff (1995) along with 
our preferred Mv(RR) relation, the Sgr dwarf is lo- 
cated at Rgc = 15.6 kpc. This rather large range 
in Rgc suggests that perhaps Ter 7, Ter 8 and Arp 
2 arc not associated with Sgr. However, there are 
considerable uncertainties associated with determin- 
ing the Galacto-centric distances; errors in the red- 
dening, metallicity, the magnitude of the HB, and 
the uncertainty in the correct Mv(RR) relation all 
lead to uncertainty in the derived Rgc distances. In 
this regard, we note that slightly different choices for 
the input parameters, can yield Galacto-centric dis- 
tances which agree within 0.9 kpc for Ter 7, Ter 8 
and Sgr, with Arp 2 still being somewhat anomalous, 
with a Rgc value which is about 3 kpc higher than 
the other objects. While a definitive answer will only 



come from proper motion studies, it appears that Ter 
7, Ter 8, and Sgr are associated. Although the evi- 
dence for Arp 2 being associated with Sgr is not as 
strong, it still remains a possibility, which should be 
investigated further. Thus, two of the anomalously 
yoTing GCs were likely formed as part of Sgr, and are 
now being accreted onto our Galaxy. 

Lin & Richer (1992) and Buonanno et al. (1994) 
have suggested that Pal 12, Arp 2. Rup 106 and Ter 7 
may all have been captured by our Galaxy, and repre- 
sent later infall events. As such, they are not indica- 
tive of the early formation of the Galactic halo. This 
argument is based on the fact that these four clusters 
lie along a single great circle, which passes through 
the Magellanic Stream. A similar argument holds for 
IC4499 (Fusi Pecci, Bellazzini & Ferraro 1995). How- 
ever, even if they have been captured by the Galaxy, 
their formation occured within the halo. It is clear 
that these young clusters formed much later than the 
majority of GCs in the Galactic halo. Thus, it is true 
that these clusters were not part of the early halo for- 
mation in the Galaxy. However, whether these clus- 
ters are later accretion events or not, they are part 
of the Galactic halo, and so give us insights into its 
formation. 

4. Age, Metallicity and Galactocentric Dis- 
tance 

The question of whether or not an age-metallicity 
relationship exists in the Galactic halo is a long 
standing problem. It has long been realized that if 
Mv(RR) has a shallow slope with metallicity, then an 
age-metallicity relationship will exist, with the most 
metal-poor clusters being the oldest (e.g. Sandage 
1982; Sarajedini & King 1989). This trend is illus- 
trated in Figure 1, which plots the age as a function 
of metallicity, assuming Mv (RR) = 0.20 [Fe/H] + 0.98. 
A least squares fit to all of the data, which takes into 
account the error bars in both coordinates (Press et 
al. 1992), yields ig ^ (-4.0±0.4) [Fe/H] -f (9.8±0.7), 
with the probability of a non-zero slope being greater 
than 99.999%. If the disk clusters (47 Tuc, NGC 6352, 
6838) and the GCs which are likely associated with 
Sgr (Ter 7, Ter 8 and Arp 2) are removed from the fit 
(to leave a pure halo sample) then the least squares 
solution yields tg = (-4.4 ± 0.9) [Fe/H] + (9.3 ± 1.4), 
with the probability of a non-zero slope being greater 
than 99.999%. If the sample is further divided by 
removing the remaining young clusters or metal-rich 
clusters ([Fe/H] > -1; NGC 6171, 6652, Rup 106, 
IC4499, and Pal 12) then we find ts = (-3.4 ± 
1.0) [Fe/H] -I- (11 ± 2), with the probability of an age- 
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Table 3 
Globular Cluster Ages 





0.17 [Fe/H] 


0.20 [Fe/H] 


0.15 [Fc/H] 


0.15 [Fc/H] 


0.20 [Fe/H] 


0.20 [Fe/H] 


0.25 [Fc/H] 


0.25 [Fc/H] 


0.30 [Fe/H] 


0.30 [Fc/H] 




+0.79 


+0.98 


+0.98 


+0.725 


+1.06 


+0.82 


+1.14 


+0.915 


+1.22 


+1.01 




A ere 


A ere 


A ere 


A ere 


A ere 


A ere 


A ere 


A ere 


A ere 


A ere 


104 


12.4 ±1.1 


15.6 ± 1.3 


16.4+1.3 


11.6 + 1.0 


17.3+1.4 


12.6+1.1 


18.4+1.5 


13.6 + 1.2 


19.4+1.5 


14.8+1.2 


288 


16.5 ± 1.3 


19.4 + 1.6 


20.9 + 1.7 


15.8 + 1.3 


21.1 + 1.7 


16.2 + 1.3 


21.4 + 1.7 


16.7 + 1.3 


21.7 + 1.7 


17.3 + 1.4 


362 


14.3 ± 1.7 


17.4 + 2.0 


18.9 + 2.1 


13.5 + 1.6 


19.3 + 2.1 


14.1 + 1.7 


19.8 + 2.2 


14.7 + 1.7 


20.3 + 2.2 


15.3 + 1.8 


1261 


13.7 ± 1.1 


16.2 + 1.3 


17.4 + 1.4 


13.1 + 1.1 


17.6 + 1.4 


13.5 + 1.1 


17.9 + 1.4 


14.0 + 1.1 


18.3 + 1.5 


14.5 + 1.2 


1851 


12.0 ±0.8 


14.2 + 1.0 


15.3 + 1.1 


11.6 + 0.8 


15.5 + 1.1 


11.9 + 0.8 


15.7+1.1 


12.3 + 0.8 


16.0 + 1.1 


12.7 + 0.8 


1904 


14.5 ± 1.3 


16.9 + 1.5 


18.5+ 1.7 


14.1 + 1.3 


18.4 + 1.7 


14.2 + 1.3 


18.3 + 1.7 


14.4 + 1.3 


18.4 + 1.6 


14.7 + 1.3 


2298 


13.9 ± 1.9 


16.0 + 2.2 


17.7 + 2.5 


13.4 + 1.8 


17.4 + 2.4 


13.5 ± 1.8 


17.2 + 2.4 


13.5 + 1.8 


17.1 + 2.3 


13.7 + 1.9 


2808 


14.7 ± 1.4 


17.3 + 1.6 


18.7+1.8 


14.1 + 1.3 


18.9 + 1.8 


14.5 ± 1.4 


19.2 + 1.8 


15.0 + 1.4 


19.5 + 1.8 


15.5 + 1.4 


3201 


11.9 ± 1.3 


13.8 + 1.6 


15.0 + 1.7 


11.5 + 1.3 


15.1 + 1.7 


11.6 ± 1.3 


15.0 ± 1.7 


11.8 + 1.3 


15.1 + 1.7 


12.1 + 1.3 


4147 


15.4 ± 1.4 


17.9 + 1.6 


19.7+1.8 


14.9 + 1.4 


19.5 ± 1.8 


15.0 ± 1.4 


19.3+ 1.7 


15.1 ± 1.4 


19.1 + 1.7 


15.3 + 1.4 


4590 


13.6 ± 0.9 


15.5 + 1.0 


17.3 + 1.2 


13.2 + 0.9 


16.8 + 1.1 


13.1 + 0.9 


16.4 + 1.1 


13.0 + 0.8 


16.0 + 1.0 


12.9 + 0.8 


5024 


15.4 ± 1.4 


17.6 + 1.6 


19.6 + 1.8 


15.0 + 1.4 


19.2 + 1.7 


14.9 + 1.3 


18.7+1.7 


14.8 + 1.3 


18.4 + 1.7 


14.8 + 1.3 


5053 


15.5 ±1.3 


17.4 + 1.5 


19.7+1.7 


15.2 + 1.3 


18.9 + 1.6 


14.9 + 1.2 


18.2 + 1.6 


14.5 + 1.2 


17.4 + 1.5 


14.1 + 1.2 


5272 


14.1 ±0.8 


16.4 ±1.0 


17.9 + 1.1 


13.6 + 0.8 


17.9 + 1.1 


13.8 + 0.8 


17.8 + 1.1 


14.0 + 0.8 


17.8 + 1.1 


14.2 + 0.8 


5466 


16.5 ±1.9 


18.7 ±2.0 


21.0 + 2.4 


16.1 + 1.9 


20.3 + 2.1 


15.9 + 1.7 


19.7 + 1.9 


15.6 + 1.5 


19.1 + 1.7 


15.4 + 1.4 


5897 


15.1 ±1.8 


17.5 ±2.1 


19.2 + 2.3 


14.6 + 1.8 


19.1 + 2.3 


14.7+1.8 


19.0 + 2.3 


14.9 + 1.8 


19.0 + 2.3 


15.2 + 1.8 


5904 


14.6 ± 1.1 


17.2 ±1.3 


18.6+1.4 


14.0 + 1.0 


18.8 + 1.4 


14.4 + 1.1 


19.0 + 1.4 


14.8 + 1.1 


19.2 + 1.4 


15.3 + 1.1 


6101 


12.5 ± 1.2 


14.4 ±1.3 


15.9 + 1.5 


12.1 + 1.1 


15.7+1.5 


12.2 + 1.1 


15.5 + 1.4 


12.3 + 1.1 


15.4 + 1.4 


12.4 + 1.1 


6121 


15.5 ±1.7 


18.3 ±2.0 


19.6 + 2.1 


14.8 + 1.6 


19.9 + 2.1 


15.3 + 1.6 


20.2 + 2.1 


15.8 + 1.7 


20.6 + 2.2 


16.4 + 1.7 


6171 


15.7± 1.9 


18.8 ±2.3 


19.9 + 2.4 


15.0 + 1.8 


20.5 + 2.5 


15.7+1.9 


21.2 + 2.5 


16.6 + 2.0 


21.9 + 2.6 


17.4 + 2.1 


6205 


16.1 ± 2.2 


18.8 + 2.6 


20.5 + 2.8 


15.5 + 2.1 


20.5 + 2.8 


15.8 ± 2.2 


20.4 + 2.8 


16.0 + 2.2 


20.4 + 2.8 


16.3 ± 2.2 


6218 


12.0 ± 1.1 


14.1 + 1.3 


15.2 + 1.4 


11.5 + 1.0 


15.4 + 1.4 


11.9 + 1.1 


15.6 + 1.4 


12.2+ 1.1 


15.9 + 1.5 


12.7 + 1.1 


6254 


17.4 ± 1.7 


20.4 + 2.0 


22.2 + 2.2 


16.8 + 1.7 


22.2 + 2.1 


17.1 + 1.7 


22.2 + 2.1 


17.4+ 1.7 


22.3 + 2.1 


17.8 + 1.7 


6341 


19.4 ± 1.5 


22.1 + 1.7 


24.8 + 1.9 


19.0 + 1.4 


24.0 + 1.8 


18.7 + 1.4 


23.3 + 1.7 


18.4+ 1.4 


22.6 + 1.7 


18.2 + 1.4 


6352 


12.9 ± 1.1 


16.3 + 1.3 


16.8 + 1.3 


12.0 + 1.0 


17.9 + 1.4 


13.2 + 1.1 


19.3 + 1.5 


14.5 + 1.2 


20.6 + 1.5 


15.9 + 1.3 


6397 


18.4 ± 1.7 


21.2 + 1.9 


23.4 + 2.2 


17.8 + 1.7 


23.0 + 2.1 


17.8 + 1.6 


22.7 + 2.0 


17.8 + 1.6 


22.4 + 2.0 


18.0 + 1.6 


6535 


16.3 ± 2.1 


18.9 + 2.4 


20.8 + 2.6 


15.8 + 2.0 


20.6 + 2.6 


15.9 ± 2.0 


20.4 + 2.5 


16.0 + 2.0 


20.4 ± 2.5 


16.3 + 2.0 


6584 


12.7 ± 1.2 


14.9 + 1.4 


16.2 + 1.6 


12.3 + 1.2 


16.2 + 1.5 


12.5 + 1.2 


16.3 + 1.5 


12.8 + 1.2 


16.4 + 1.5 


13.1 + 1.2 


6652 


9.1 ± 1.3 


11.4 + 1.7 


12.2 + 1.8 


8.5 + 1.2 


12.8 + 1.8 


9.1 + 1.3 


13.5 + 1.9 


9.8 + 1.4 


14.1 + 1.9 


10.5 + 1.5 


6752 


17.7 ± 1.9 


20.7 + 2.2 


22.5 + 2.3 


17.0 + 1.8 


22.5 + 2.3 


17.4 + 1.8 


22.6 + 2.3 


17.7 + 1.8 


22.8 + 2.3 


18.2 + 1.9 


6809 


16.5 ± 1.2 


19.1 + 1.3 


21.0 + 1.5 


16.0 + 1.2 


20.8 + 1.4 


16.1 + 1.1 


20.5 + 1.4 


16.1 + 1.1 


20.4 + 1.3 


16.3 + 1.1 


6838 


11.3 ±0.9 


14.3 + 1.1 


14.8 + 1.1 


10.5 + 0.8 


15.8 + 1.1 


11.5 + 0.9 


16.9 + 1.2 


12.5 + 0.9 


18.0 + 1.3 


13.7 + 1.0 


7006 


14.0 ± 1.1 


16.4 + 1.3 


17.9 + 1.4 


13.5 + 1.1 


17.9 + 1.4 


13.8 + 1.1 


17.9 + 1.4 


14.0 + 1.1 


18.0 + 1.4 


14.3 + 1.1 


7078 


17.1 ±1.7 


19.5 ±2.0 


21.9 + 2.2 


16.7 + 1.7 


21.2 + 2.2 


16.5 + 1.7 


20.7 + 2.1 


16.3 + 1.7 


20.2 + 2.1 


16.2 + 1.7 


7099 


16.9 ± 1.6 


19.3 ±1.8 


21.6 + 2.0 


16.5 + 1.5 


21.0 + 1.9 


16.3 + 1.5 


20.4 + 1.8 


16.1 + 1.5 


19.9 + 1.8 


16.0 + 1.4 


7492 


17.6 ± 1.9 


20.3 ±2.0 


22.4 + 2.4 


17.0 + 1.9 


22.1 + 2.2 


17.1 + 1.7 


21.9 + 2.1 


17.2 + 1.6 


21.7 + 2.0 


17.4 + 1.6 


Ter 7 


7.2 ±0.5 


8.7 ±0.8 


9.0 + 0.8 


6.8 + 0.4 


9.7 + 0.9 


7.3 + 0.5 


10.6 + 1.0 


7.9 + 0.7 


11.4 + 1.1 


8.7 + 0.8 


Ter 8 


16.9 ± 1.5 


19.4 ±1.7 


21.5 + 1.9 


16.4 + 1.5 


21.1 + 1.8 


16.4 + 1.4 


20.7+1.8 


16.3 + 1.4 


20.4 + 1.8 


16.3 + 1.4 


Rup 106 


13.2 ±0.8 


15.7 ±0.9 


17.4+1.0 


12.7 + 0.8 


17.4+1.0 


12.9 + 0.8 


17.3 + 1.0 


13.0 + 0.8 


17.2 + 1.0 


13.2 + 0.8 


Pal 5 


13.4 ± 1.4 


15.7 + 1.6 


17.0 + 1.9 


12.9 + 1.4 


17.1 ± 1.8 


13.2 + 1.4 


17.2 + 1.7 


13.5 + 1.3 


17.4 + 1.6 


13.9 + 1.3 


Pal 12 


9.4 ± 1.4 


11.7 + 1.7 


12.6 + 1.9 


8.8 + 1.4 


13.1 ± 1.8 


9.3 + 1.4 


13.5 + 1.8 


9.8 + 1.4 


13.9 ± 1.8 


10.4 + 1.4 


IC4499 


10.6 ± 1.1 


12.2 + 1.2 


13.4 + 1.4 


10.3 + 1.0 


13.3 + 1.3 


10.4 + 1.0 


13.1 + 1.3 


10.4 + 1.0 


13.1 ± 1.3 


10.6 + 1.0 


Arp 2 


12.3 ±0.8 


14.2 + 1.0 


15.6 + 1.1 


11.9 + 0.8 


15.5 + 1.1 


12.0 + 0.8 


15.4+1.0 


12.1 + 0.8 


15.4 + 1.0 


12.4 + 0.8 
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metallicity relation existing being greater than 99.7%. 
Note that the above results are contingent upon a 
constant [a/Fe] below [Fe/H] = —1.0, as this was as- 
sumed in the isochrones. However, we note that to re- 
move the age-metallicity trend observed would require 
that [a/Fe] increase by 0.8 dex, from [Fe/H] = —2.2 
to [Fe/H] = —1.2, the [Fe/H] range spanned by the 
majority of the GCs in the sample. However, such a 
dramatic change in [a/Fe] is not consistent with the 
halo star observations, as discussed in §2.1. 

To further explore the age-metallicity question, GC 
ages have been determined using our best estimate for 
the Mv(RR) zero-point, and slopes ranging from 0.15 
to 0.30, in steps of 0.002. The resulting ages were then 
analyzed using the least-squares fit as above, in order 
to examine how the slope of Mv(RR) with metallic- 
ity affects the age-metallicity relationship. The re- 
sults are shown in Figure 2, which plots the proba- 
bility of an age-metallicity relation as a function of 
the Mv(RR) slope with metallicity. If all clusters are 
included in the analysis, then an age-metallicity rela- 
tion exists at the greater than 99.9% confidence level 
for all values of the slope tested (from 0.15 to 0.30). 
If the disk clusters and the GCs associated with Sgr 
are removed from the fit, leaving a pure halo sam- 
ple, then an age-metallicity relationship exists at a 
greater than 97.6% confidence level for all slopes less 
than or equal to 0.3. If the young and/or metal-rich 
clusters are excluded, then an age-metallicity rela- 
tionship exists at the 2 <j (95%) level, provided that 
the slope of Mv(RR) with metallicity is less than 
0.259. As most evidence favors low values for the 
Mv(RR) slope (see also §6), Figure 2 demonstrates 
that an age-metallicity relationship exists in the halo 
of our Galaxy. 

The relationship between age and Galactocentric 
distance is shown in Figure 3, which plots age (as- 
suming Mv(RR) = 0.20 [Fe/H] -|- 0.98) as a function 
of Rgc- The sample has been classified into old halo, 
younger halo and disk clusters based on metallicity, 
kinematics and HB morphology (Zinn 1993). This 
will be discussed in more detail in section 5. Here 
we note that no clear relationship between age and 
Galactocentric distance exists, though there is a sug- 
gestion that GCs become younger as one goes to large 
Galactocentric distances. A least-squares fit to the 
data yields tg = (-0.06±0.03) Rgc + (16.5±0.7) with 
non-zero slope being significant only at the 94.5% 
confidence level. If the apogalactica distances given 
in Table 2 are substituted for Rqc where available, 
then the significance of the non-zero slope drops below 
50%. Using ages and distances derived from the other 
HB relationships given in Table 1 results in similar 



Fig. 1. — Age as a function of metallicity for 43 GCs, 
assuming Mv(RR) = 0.20 [Fe/H] 0.98. 



Fig. 2. — The probability of an age-metallicity rela- 
tion is plotted as a function of the slope with metal- 
licity of Mv(RR). If all clusters are included (not 
plotted), then the probability is greater than 0.999 
for all values of the slope. The halo cluster sample 
does not include 47 Tuc, NGC 6352, NGC 6838, Ter 
7, Ter 8 and Arp 2, while the old halo sample excludes 
NGC 6171, 6652, Rup 106, IC4499, and Pal 12 in ad- 
dition to the above clusters. Provided that the slope 
of Mv(RR) with metallicity is less than 0.26, an age- 
metallicity relationship exists in the halo, regardless 
of which sample is used to define the halo. 
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Fig. 3. — Age as a function of Galactoccntric distance 
(Rgc) for 43 GCs, assuming Mv(RR) = 0.20 [Fc/H] + 
0.98. The sample has been classified into old halo, 
younger halo and disk clusters based on metallicity, 
kinematics and HB morphology (Zinn 1993). Error 
bars have not been plotted for clarity, but are typi- 
cally ±1.6 Gyr. 

plots. Hence, the present data provide no compelling 
evidence for an age-Galactocentric distance relation- 
ship. 

5. Age Range 

In Table 3 there are a wide range of GC ages for 
a given Mv(RR) relation. Clearly, some of this is 
due to the relatively large errors (of order ±1.6 Gyr) 
in the individual age determinations. To quantify 
how much of the age range is due to the errors, and 
whether an intrinsic age range exists within the GC 
system, the following statistical test was performed: 
an 'expected' distribution for no intrinsic age range 
was constructed by randomly generating 10,000 ages 
using a Gaussian distribution, with a mean given by 
the mean age of the entire sample, and the signia 
(i.e. standard deviation) given by the error in an in- 
dividual age determination. This is repeated for all 
clusters in the sample, so that the expected distri- 
bution contains 43 x 10, 000 = 430, 000 ages. This 
expected distribution is then compared to the actual 
age distribution, using the F-test (Press et al. 1992), 
which determines if the two distributions have the 
same variance. If there is less than a 5% chance that 
the two distributions have the same variance, then 
we conclude that an age range exists. The size of 
the age range is inferred by the standard method, 



Fig. 4. — Histogram of ages for our preferred 
Mv(RR) relation. The solid line is the expected his- 
togram of ages given the errors in the individual ages, 
and assTiming no intrinsic age range. It has been nor- 
malized to the total number of clusters in our sample 
(43). It is clearly not a good fit to the data, and the 
F-test rejects the hypothesis that the two distribu- 
tions have the same variance at a very high confidence 
level. The dotted line shows the best fitting Gaussian 
distribution, which includes an intrinsic age range of 
frange = 2.3 Gyr, in addition to the scatter induced 
by the error in the individual ages. 



(Grange = (^l^^ - (^^pected' ^^ere a^ange is the sigma 
of the true age range, Uohs is the sigma of the actual 
data, and (Toxpoctod is the sigma of our expected dis- 
tribution, given the input errors in our ages. Tests 
have been performed which indicated that the typical 
error in our inferred Urange is ±0.1 Gyr for a given 
Mv(RR) relation. 

When this analysis is performed on all 43 clus- 
ters, the F-test rejects the hypothesis of no intrin- 
sic age range, at a very high confidence level for all 
Mv(RR) relations used in this paper. Figure 4 plots 
the actual age histogram for our preferred Mv(RR) re- 
lation, along with the expected histogram if there was 
no age range, and the best fitting histogram, which 
includes an intrinsic range of ages, with Urange = 2.3 
Gyr. If the age spread is defined to be the age 
range which includes 95% of the clusters, then the 
age spread inferred is 4 x 2.3 = 9.2 Gyr. The com- 
plete results of this analysis are presented in detail in 
Table 4, for all of the Mv(RR) relations. 

In looking at Figure 4, it is clear that the 4 very 
young clusters (Ter 7, Pal 12, IC 4499 and NGC 



10 



TABLE 4 
AGE RANGE 



Probability of (Grange Age 
Mv(RRLyr) NO Age Range (Gyr) Spread 

All Clusters (N=43) 



0.17[Fc/H] - 


F0.79 


2.1 X 10" 


-12 


2.2 


8.7 


0.20[Fc/H] - 


F0.98 


1.7 X 10" 


-9 


2.3 


9.2 


0.15[Fc/H] - 


F0.98 


1.8 X IQ- 


-11 


2.7 


10.8 


0.15[Fc/H] - 


1-0.72 


6.9 X 10- 


-14 


2.2 


8.9 


0.20[Fc/H] - 


1- 1.06 


2.6 X 10- 


-9 


2.4 


9.8 


0.20[Fc/H] - 


1-0.82 


4.2 X 10- 


-11 


2.0 


8.1 


0.25[Fe/H] - 


1- 1.14 


2.2 X 10" 


-7 


2.2 


8.8 


0.25[Fe/H] - 


1-0.91 


5.5 X 10" 


-9 


1.9 


7.5 


0.30[Fe/H] - 


1- 1.22 


2.3 X 10" 


-6 


2.1 


8.3 


0.30[Fe/H] + 1.01 


2.8 X 10- 


-7 


1.8 


7.0 



Excluding Young Clusters (N=39) 



0.17[Fc/H] - 


1-0.79 


7.1 X 10-" 


1.4 


5.7 


0.20[Fc/H] - 


F0.98 


9.0 X 10-3 


1.4 


5.6 


0.15[Fc/H] - 


F0.98 


1.2 X 10-3 


1.8 


7.1 


0.15[Fc/H] - 


F0.72 


1.2 X 10-" 


1.5 


6.0 


0.20[Fc/H] - 


h 1.06 


1.4 X 10-2 


1.5 


5.9 


0.20[Fc/H] - 


hO.82 


3.3 X 10-3 


1.3 


5.1 


0.25[Fe/H] - 


h 1.14 


4.1 X 10-2 


1.3 


5.2 


0.25[Fc/H] - 


F0.91 


2.0 X 10-2 


1.1 


4.4 


0.30[Fc/H] - 


1- 1.22 


4.7 X 10-2 


1.3 


5.1 


0.30[Fe/H] - 


1- 1.01 


4.3 X 10-2 


1.0 


4.1 



Clusters with -1.8 < [Fe/H] < -1.3 (N=21) 



0.17[Fc/H] - 


1-0.79 


1.9 X 10-2 


1.3 


5.3 


0.20[Fc/H] - 


hO.98 


1.0 X 10-2 


1.6 


6.4 


0.15[Fc/H] - 


F0.98 


1.9 X 10-2 


1.7 


6.8 


0.15[Fc/H] - 


F0.72 


1.6 X 10-2 


1.3 


5.1 


0.20[Fc/H] - 


^ 1.06 


1.3 X 10-2 


1.7 


6.9 


0.20[Fc/H] - 


1-0.82 


1.9 X 10-2 


1.3 


5.2 


0.25[Fc/H] - 


1- 1.14 


1.0 X 10-2 


1.8 


7.0 


0.25[Fc/H] - 


f-0.91 


1.3 X 10-2 


1.4 


5.4 


0.30[Fe/H] - 


1- 1.22 


7.7 X 10-3 


1.8 


7.2 


0.30[Fe/H] - 


\- 1.01 


6.0 X 10-3 


1.4 


5.8 



6652) are somewhat anomalous, and responsible for a 
good part of the very large inferred age range. Thus, 
the above analysis was repeated excluding the 4 very 
young clusters. This gives a reasonable estimate of 
the true age range among the bulk of the Galac- 
tic GCs. Even with this restricted sample, an in- 
trinsic age range exists at the greater than the 95% 
confidence for all Mv(RR) relations (Table 4). The 
size of the age range is reduced, and varies between 
4.1 7.1 Gyr, depending on the choice of Mv(RR). 
The sensitivity of these results to the slope of the 
Mv (RR) relation with metallicity may be reduced by 
considering only those GCs in the restricted range 
-1.8 < [Fe/H] < -1.3. There are 21 GCs in this 
group (the young clusters Ter 7, Pal 12 and NGC 
6652 are not included), and the results are quite simi- 



lar to those obtained with the sample which excludes 
the very young clusters. As shown in Table 4, an 
intrinsic age range exists regardless of the choice of 
Mv(RR). The age spread is 5.1 — 7.2, depending on 
the choice of Mv(RR). Thus, we may conclude that 
a real age spread of ^ 5 Gyr exists among the bulk 
of the GCs, with several clusters (~ 10%) which are 
considerably younger. 

6. The Second Parameter Problem 

The morphology of the HB (ie. the relative num- 
ber of red, blue and RR Lyr stars on the HB) is pri- 
marily governed by the metallicity of the cluster. As 
such, [Fe/H] is the first parameter which controls HB 
morphology. However, it has long been known that 
two clusters with similar [Fe/H] values, can have con- 
siderably different HB morphologies. NGC 288 and 
362, and M13 and M3 are classic examples of GCs 
which demonstrate that some other parameter besides 
[Fe/H] is important in determining the morphology 
of the HB. Searle & Zinn (1978) demonstrated that 
the second parameter is correlated with Galactoccn- 
tric distance; there is a tight relationship between 
[Fe/H] and HB type in the inner halo (Rgc < 8kpc), 
while the effects of the second parameter are most 
pronounced in the outer halo. The quest to deter- 
mine the nature of the second parameter which gov- 
erns HB morphology has been a longstanding one in 
astronomy. There are numerous possibilities for the 
second parameter (age, oxygen abundance, core ro- 
tation, mass loss on the RGB, etc.). Given that the 
previous section has demonstrated that a large in- 
trinsic age range exists among the GCs, we will focus 
here on examining the hypothesis that age is the sec- 
ond parameter (Searle & Zinn 1978; Lee, Demarque 
& Zinn 1994). 

On the assumption that age is the second param- 
eter, Zinn (1993) has divided the halo GCs into two 
groups, the Old Halo (OH) and Younger Halo (YH; 
these groupings arc given in Table 2). GCs were 
deemed to be younger if their HB types were 0.4 red- 
der (using the (B — R)/(B-|-V-|-R) index, see footnote 
3 in §2.3) than the typical inner halo cluster at their 
metallicity. There are 25 OH clusters in our sample, 
and 15 young halo clusters. Of the clusters which are 
clearly young in our sample. IC4499, Rup 106, Pal 12 
and Ter 7 are all part of the YH grouping. Only Arp 
2 is incorrectly classified as a OH cluster. In addition, 
NGC 6652 has a a young AV^g age, even though it 
is classified as a OH cluster. This suggests that age is 
the dominant second parameter. Indeed, the YH clus- 
ters do tend to have lower ages than the OH clusters. 
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as shown in Figure 3. However, incorrect classifica- 
tion of Arp 2 suggests that a third parameter affects 
the HB type of some clusters. Note that NGC 6652 
([Fe/H] = -0.89 and HB type = -1.00) lies at the 
boundary of the OH and YH clusters, so whether it 
belongs to the OH or YH group is uncertain. 

The errors in our AV^g ages can be rather large, 
thus, for the bulk of the GCs, it is difficult to say with 
certainty that one particular GC is younger than an- 
other. This difficulty may be overcome by determin- 
ing the weighted mean age of the OH and YH groups. 
The results of this calculation are shown in Table 5. 
If all of the halo clusters are included in the sam- 
ple, then the YH is 2 - 4 Gyr younger than the OH 
group, and the difference in age is significant at the 
4.8 — 7.5 a level. Perhaps more importantly, if the GCs 
are randomly sorted into two groups of the same size 
as the YH and OH groups, age differences as large as 
those found between the YH and OH groups only oc- 
cur 0.5% of the time. Given the spread in ages found 
in the previous section, this latter test ensures that 
the differences in the mean ages of the two groups is 
not just a coincidence. If the Sgr clusters (Ter 7, Ter 
8 and Arp 2) are removed, then an age difference of 
2.5 Gyr is found significant at the 4.3 — 4.6 a level. 
The chance of such a large age difference occurring 
in random subgroups is less than 0.5%. If, in addi- 
tion to the Sgr clusters, all of the young clusters are 
excluded from the sample (IC4499, Pal 12, Rup 106, 
NGC 6652), then the difference in age drops to 1.6 - 
2.3 Gyr (at the 3.2 — 3.5 cr level), depending on the 
choice of the Mv(RR) relation. The chance of such a 
large age difference occurring in random subgroups is 
less than 2.0%. Thus, we see that even when the ob- 
viously young clusters are removed from the sample, 
there is still a significant difference in the mean age of 
the OH and YH groups. Together, these results are 
consistent with the hypothesis that age is the second 
parameter, and that a typical second parameter clus- 
ter is about 2-3 Gyr younger than the clusters which 
possess bluer HBs at similar metallicities. 

Although our results are consistent with age being 
the second parameter, they cannot entirely rule out 
other phenomena. For example, we have assumed 
that the helium abundance is the same for all clus- 
ters of a given metallicity. If two clusters of the 
same metallicity have different helium abundances, 
it is possible to mimic the effects of a more youth- 
ful AV^B age, while reddening the HB morphology^. 
Thus, our results cannot conclusively prove that age 

^This example was pointed out to us by the referee, Bruce 
Carney. 



TABLE 5 

MEAN AGES OF THE OLD & YOUNGER HALO 









Younger 










Old Halo 


Halo 


AAgo 


AAge 


Mv(RRLyr) 


Age (Gyr) 


Age (Gyr) 


(Gyr) 


a 


All Halo Clusters 


0.17[Fe/H] - 


f 0.79 


14.9 ± 0.3 


11.6 ± 0.3 


3.3 ± 0.4 


7.5 


0.20[Fe/H] - 


f 0.98 


17.4 ± 0.4 


14.2 ± 0.4 


3.2 ± 0.5 


6.2 


0.15[Fe/H] - 


h 0.98 


19.1 ± 0.4 


15.3 ± 0.4 


3.8 ± 0.6 


6.6 


0.15[Fe/H] - 


h 0.72 


14.4 ± 0.3 


10.8 ± 0.3 


3.6 ± 0.4 


8.3 


0.20[Fe/H] - 


f 1.06 


19.0 ± 0.4 


15.7 ± 0.4 


3.3 ± 0.6 


5.9 


0.20[Fe/H] - 


h 0.82 


14.6 ± 0.3 


11.6 ± 0.3 


3.0 ± 0.4 


7.0 


0.25[Fe/H] - 


f 1.14 


18.9 ±0.4 


15.9 ± 0.4 


3.0 ±0.6 


5.4 


0.25[Fe/H] - 


f 0.91 


14.8 ±0.3 


12.2 ± 0.3 


2.6 ±0.4 


6.0 


n "^nrPi^ /rri 

u.ou[re/ nj 


f 1.22 


18.8 ±0.4 


16.2 ± 0.4 


2.6 ±0.5 


ZL ^ 


U.oU[re/rlJ - 


f 1.01 


15.0 ±0.3 


12.7 ±0.3 


2.3 ±0.4 


O.v) 






Excluding 


Sgr Clusters 






0.17[Fc/H] - 


f 0.79 


15.2 ±0.3 


13.0 ±0.3 


2.2 ±0.5 


4.6 


0.20[Fe/H] - 


f 0.98 


17.8 ±0.4 


15.3 ± 0.4 


2.5 ±0.5 


4.5 


0.15[Fe/H] - 


f 0.98 


19.5 ±0.4 


16.7 ±0.4 


2.7 ±0.6 


4.5 


0.15[Fe/H] - 


f 0.72 


14.6 ±0.3 


12.5 ± 0.3 


2.1 ±0.5 


4.6 


0.20[Fe/H] - 


f 1.06 


19.4 ±0.4 


16.7 ±0.4 


2.7 ±0.6 


4.5 


0.20[Fe/H] - 


f 0.82 


14.9 ±0.3 


12.8 ± 0.3 


2.1 ±0.5 


4.6 


0.25[Fe/H] - 


f 1.14 


19.3 ±0.4 


16.7 ±0.4 


2.6 ±0.6 


4.4 


0.25[Fe/H] - 


f 0.91 


15.1 ±0.3 


13.0 ±0.3 


2.1 ±0.5 


4.6 


0.30[Fe/H] - 


f 1.22 


19.2 ±0.4 


16.8 ± 0.4 


2.5 ±0.6 


4.3 


0.30[Fe/H] - 


f 1.01 


15.3 ±0.3 


13.2 ± 0.3 


2.1 ±0.5 


4.5 




Excluding Sgr & 


Very Young Clusters 




0.17[Fe/H] - 


f 0.79 


15.2 ±0.3 


13.5 ± 0.4 


1.7±0.5 


3.4 


0.20[Fe/H] - 


f 0.98 


17.8 ±0.4 


15.7 ±0.4 


2.0 ±0.6 


3.4 


0.15[Fe/H] - 


f 0.98 


19.5 ±0.4 


17.1 ± 0.5 


2.3 ±0.7 


3.5 


0.15[Fe/H] - 


f 0.72 


14.6 ±0.3 


13.0 ±0.4 


1.7±0.5 


3.4 


0.20[Fe/H] - 


f 1.06 


19.4 ±0.4 


17.2 ± 0.5 


2.2 ±0.6 


3.5 


0.20[Fe/H] - 


f 0.82 


14.9 ±0.3 


13.2 ± 0.4 


1.7±0.5 


3.4 


0.25[Fe/H] - 


f 1.14 


19.3 ±0.4 


17.2 ± 0.5 


2.1 ±0.6 


3.3 


0.25[Fe/H] - 


f 0.91 


15.1 ±0.3 


13.4 ±0.4 


1.7±0.5 


3.3 


0.30[Fe/H] - 


f 1.22 


19.2 ±0.4 


17.3 ± 0.5 


2.0 ±0.6 


3.2 


0.30[Fe/H] - 


f 1.01 


15.3 ±0.3 


13.7 ±0.4 


1.6 ±0.5 


3.2 



is the second parameter, only that our ages are con- 
sistent with age being the second parameter. How- 
ever, we note that Lee et al. (1994) have extensively 
discussed arguments against parameters besides age 
being responsible for the second parameter. They 
found problems with every candidate second parame- 
ter except age. For example, variations in the helium 
abundance are ruled out by constraints set by RR Lyr 
periods. 

As mentioned above, the second parameter is cor- 
related with Galactocentric distance, but, as demon- 
strated by Figure 3, and discussed in the previous 
section, age is not correlated with Galactocentric dis- 
tance in our data. This would appear to contradict 
the conclusion that age is the second parameter. This 
seeming contradiction may be resolved by a few fac- 
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tors: (1) the relatively small number of GCs in our 
sample coupled with the relatively large age errors 
(Lee et al. 1994 have 83 GCs in their HB sample); (2) 
some of our clusters arc far from their apogalacticon, 
which tends to reduce the size of radial gradients; and 
(3) in the region Rgc =8 — 40 kpc, there is consider- 
able scatter in the HB type-[Fe/H] correlation, with 
some clusters in this region having a similar HB type 
to clusters in the inner halo (Rgc = 8 kpc) . This is 
illustrated by the fact that beyond 8 kpc, our sample 
contains 13 OH clusters, and 14 YH clusters. This 
suggests that while we should find a greater range 
of ages in the outer halo sample, the oldest clusters 
in the outer halo will have a similar age to the old- 
est clusters in the inner halo. Hence, one does not 
expect to find a strong age-Galactocentric distance 
relationship even if age is the second parameter. In- 
stead, there should be a greater age range in the outer 
halo, as compared to the inner halo. To test this hy- 
pothesis, the age range calculations discussed in the 
previous section were applied to the inner halo and 
outer halo sample. When all clusters were included, 
the probability of an age range existing was much 
higher in the outer halo sample, as opposed to the in- 
ner halo sample. However, the age range of the inner 
halo sample was nearly the same as the outer halo 
sample (9.6 vs. 8.5 Gyr for our preferred Mv (RR) re- 
lation). The determination of the age range of the 
inner clusters is rather uncertain, as there are only 
13 clusters in this group, of which NGC 6652 is an 
obvious outlier. If it is removed from the inner halo 
sample, then there is no evidence for an age spread 
among the inner halo clusters. Hopefully, the ques- 
tion of an age-Galactocentric distance relation, and 
whether there is a difference in the range of ages found 
in the inner and outer halo will be resolved by more, 
high quality data. 

7. Discussion 

In analyzing our AV^g ages of 43 GCs (presented 
in Table 3), the following conclusions were drawn: 

(1) if the slope of Mv(RR) with metallicity is less 
than 0.26, then an age-metallicity relationship exists, 
with the most metal-poor clusters being the oldest; 

(2) our data set does not contain strong evidence for 
an age-Roc distance relationship; (3) independent of 
the choice of Mv (RR) there is strong evidence for an 
age spread of 5 Gyr among the bulk of the GCs; (4) 
about 10% of the GCs are substantially younger than 
the majority and including them in the total sample 
increases the age range to about 9 Gyr; and (5) the 
mean age of the red-HB, second parameter clusters 
is 2 -3 Gyr younger than normal clusters, which is 



consistent with age being the second parameter. It 
would appear that conclusions (2) and (5) contradict 
each other, since the second parameter is correlated 
with Galactocentric distance (Searle & Zinn 1978; Lee 
et al. 1994). This contradiction may be resolved by 
noting that the errors in the age determinations are 
rather large (^ 10%); there is considerable scatter in 
the HB-Rgc relation; and there is weak evidence for 
an age-RGc distance relationship (at the 94.5% con- 
fidence limit) in our data set. 

In addition to the above, we note that the GCs Ter 
7, Ter 8 and Arp 2 appear to be associated with the 
recently discovered Sgr dwarf spheroidal galaxy. This 
is based on the fact that the above objects are all 
located at similar distances and in the same region of 
the sky as Sgr (Ibata et al. 1994). In addition, these 
GCs have radial radial velocities similar to that of Sgr 
(Da Costa & Armandroff 1995). It appears that the 
Galaxy is in the process of accreting the Sgr dwarf 
and its accompanying GCs, two of which (Ter 7 and 
Arp 2) are anomalously young compared to the bulk 
of the GCs in the Galaxy. 

The above conclusions strengthen the original pro- 
posal by Scarle & Zinn (1978) that the outer halo of 
the Galaxy formed in a slow, rather chaotic collapse, 
with the Galaxy accreting material over several Gyr. 
The Sgr dwarf, and associated GCs is an example of a 
large gas fragment which collapsed, and self-enriched 
and is now being accreted by the Galaxy. However, 
not all of the outer halo (Rgc > 8 kpc) formed via 
later accretion. A significant fraction (50%) of the 
outer halo clusters in our sample do not have a strong 
second parameter effect, and so are part of Zinn's 
(1993) Old Halo group. These old, outer halo ob- 
jects were formed during the prompt collapse of the 
proto-Galactic cloud, though they still may have been 
accreted at a later time. As time passed, more metal- 
rich GCs formed, and were accreted into the outer 
halo, leading to the observed age-metallicity relation- 
ship. 

In contrast to our earlier work on the ages of 32 
GCs (Chaboyer et al. 1992), the present GC sam- 
ple contains no evidence that the inner halo formed 
in a rapid collapse (see Figure 3). However, if NGC 
6652 is removed from the sample, then there is no ev- 
idence for an intrinsic age range among the inner halo 
clusters, which would suggest that the inner halo did 
indeed form in a rapid collapse. In contrast, the con- 
clusion that a large age range exists among the outer 
halo GCs is a robust statement. A more definitive 
answer on whether or not the inner halo formed in a 
rapid collapse requires more high quality data of inner 
halo clusters. 
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A. Estimating the Gaussian Error in AV^b 

The errors in AV^g given in Table 2 are those 
quoted by the observers, who rarely specify how they 
have determined the error bar. For the statistical 
analysis which comprises the bulk of the paper, Gaus- 
sian 1-a error bars are required. Perhaps the best 
way to estimate the Gaussian error in a measurement 
is to analyze repeated observations. In this vein, we 
have searched the literature for independent measure- 
ments of AVj^B for which an error is also included. 
The results of our literature search are presented in 
Table 6, which gives the cluster name, AV^g value 
and the reference. Using the AV^^ values and er- 
rors, we have calculated and tabulated the quantity 
5 = (AYa - AVb)/(e2 + elfl"^, where AVa is the 
measured AV^g with its error (ca) as reported by 
observer a, and AV;, and Cf, are the same quantities 
reported by observer h. Essentially, 5 is simply the 
difference in the AV^g observations, normalized by 
the quoted errors. If the observers are quoting Gaus- 
sian l-cr error bars, then 5 should have a Gaussian 
distribution, with a = 1. There are 17 measurements 
of 5 in Table 6; for this sample size one would ex- 
pect 5 values of 5 in excess of 1 for Gaussian errors. 
However, this occurs only twice. This suggests that 
the reported errors are an overestimate of the Gaus- 
sian one-(T error bars. Indeed, the F-tcst (Press et al. 
1992) finds that there is only a 2% chance that 5 has 
a standard deviation of 1.0. The quantity 5 has an 
actual standard deviation of 0.61. Thus, multiplying 
the quoted errors by 0.61 yields the best estimate the 
Gaussian l-cr error in AV^g. 

Another estimate for the Gaussian l-cr error in 
AVj^B may be obtained by computing the standard 
deviation (using the small sample formulae of Keep- 
ing 1962) for each set of AV^g values given in Table 
6. This standard deviation is then compared to the 
mean AVj^ error quoted by the observers in Table 
6. Dividing the standard deviation by the mean er- 
ror and taking the average of this ratio yields 0.55. 
This is the amount by which one should multiply the 
quoted AV^g errors in order to obtain a Gaussian 



TABLE 6 

INDEPENDANT AV^g OBSERVATIONS 



Cluster 




<5 


Reference 


NGC 104 


3.61 ± 0.10 
3.81 ±0.18 
3.76 ±0.10 


-0.971 
-1.060 


Table 2 
CSD 

SK 


NGC 288 


3.73 ±0.12 
3.70 ±0.14 
3.62 ±0.10 


0.163 
0.704 


Table 2 
SK 

Bergbush 1993 
Pound et al. 1987 


NGC 1261 


3.57 ±0.12 
3.30 ±0.14 


0.379 


Table 2 

Alcaino et al. 1992b 


NGC 1851 


3.45 ± 0.10 
3.34 ±0.10 


0.778 


Table 2 
CSD 


NGC 3201 


3.45 ± 0.21 
3.44 ±0.12 


0.041 


Brewer et al. 1993 
Alcaino et al. 1989 
Cacciari 1984 


NGC 4590 


3.42 ± 0.10 
3.49 ±0.12 
3.42 ±0.10 


-0.448 
0.000 


Table 2 
CSD 

Alcaino et al. 1990 
Harris 1975 


NGC 5897 


3.60 ±0.18 
3.52 ±0.14 


0.351 


Table 2 

Ferraro et al. 1992b 


NGC 6121 


3.68 ±0.16 
3.52 ±0.10 
3.45 ±0.13 


0.848 
1.116 


Table 2 
SK 

Kanatas et al. 1995 


NGC 6171 


3.75 ± 0.18 
3.70 ±0.11 


0.237 


Table 2 

Ferraro et al. 1991 


NGC 6752 


3.77 ±0.16 
3.72 ±0.14 


0.235 


Table 2 
SK 


NGC 6809 


3.66 ±0.10 
3.54 ±0.14 


0.697 


Table 2 

Alcaino et al. 1992a 


NGC 7492 


3.72 ±0.14 
3.61 ±0.10 


0.778 


Table 2 
SK 


Rup 106 


3.32 ±0.07 
3.27 ±0.12 


0.360 


Table 2 
CSD 



1-a error. This value is quite similar to the 0.61 ob- 
tained above. To be conservative, the value of 0.61 
will be used. 
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Fig. 2. — The probability of an age-metallicity relation is plotted as a function of the slope with metallicity of 
Mv(RR). If all clusters are included (not plotted), then the probability is greater than 0.999 for all values of the 
slope. The halo cluster sample does not include 47 Tuc, NGC 6352, NGC 6838, Ter 7, Ter 8 and Arp 2, while 
the old halo sample excludes NGC 6171, 6652, Rup 106, IC4499, and Pal 12 in addition to the above clusters. 
Provided that the slope of Mv(RR) with metallicity is less than 0.26, an age-metallicity relationship exists in the 
halo, regardless of which sample is used to define the halo. 
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Fig. 3. — Age as a function of Galactocentric distance (Rqc) for 43 GCs, assuming Mv(RR) = 0.20 [Fe/H] + 0.98. 
The sample has been classified into old halo, younger halo and disk clusters based on metallicity, kinematics and 
HB morphology (Zinn 1993). Error bars have not been plotted for clarity, but are typically ±1.6 Gyr. 
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Fig. 4. — Histogram of ages for our preferred Mv(RR) relation. The solid line is the expected histogram of ages 
given the errors in the individual ages, and assuming no intrinsic age range. It has been normalized to the total 
number of clusters in our sample (43). It is clearly not a good fit to the data, and the F-test rejects the hypothesis 
that the two distributions have the same variance at a very high confidence level. The dotted line shows the best 
fitting Gaussian distribution, which includes an intrinsic age range of Urange = 2.3 Gyr, in addition to the scatter 
induced by the error in the individual ages. 
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